W.J. Metzger et al. 



1 



o 
o 

(N 
> 

o 

:^ 

00 



> 



o 

Oh 

CD 



X 



Parametrization of Bose-Einstein Correlations and Reconstruction of the Source Function in 

Hadronic Z-boson Decays using the L3 Detector* 
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Bose-Einstein correlations of pairs of identical charged pions produced in hadronic Z decays are analyzed in 
terms of various parametrizations. A good description is achieved using a Levy stable distribution in conjunction 
with a hadronization model having highly correlated configuration and momentum space, the t-model. Using 
these results, the source function is reconstructed. 

Keywords; Bose-Einstein Con'elations 



I. INTRODUCTION 

In particle and nuclear physics, intensity interferometry 
provides a direct experimental method for the determination 
of sizes, shapes and lifetimes of particle-emitting sources (for 
reviews see iH Q. H H Etl)- In particular, boson interfer- 
ometry provides a powerful tool for the investigation of the 
space-time structure of particle production processes, since 
Bose-Einstein correlations (BEC) of two identical bosons re- 
flect both geometrical and dynamical properties of the particle 
radiating source. 

Here we study BEC in hadronic Z decay. We inves- 
tigate various static parametrizations in terms of the four- 
momentum difference, Q — \/—{p\ — pi)^ and find that none 
give an adequate description of the Bose-Einstein correlation 
function. However, within the framework of models assuming 
strongly correlated coordinate and momentum space a good 
description is achieved. We then reconstruct the complete 
space-time picture of the particle emitting source in hadronic 
Z decay. 

The data used in the analysis were collected by the L3 de- 
tector 10, S [M S [Id!] at an e+e^ center-of-mass energy of 
~ 91.2 GeV. Approximately 36 million like-sign pairs of 
well-measured charged tracks of about 0.8 million hadronic Z 
decays are used [11]. 

We perform analyses on the complete sample as well as 
on two- and three-jet samples. The latter are found using 
calorimeter clusters with the Durham jet algorithm lfT2i [Tsl 
with a jet resolution parameter ycut = 0.006. To determine 
the thrust axis of the event we also use calorimeter clusters. 



II. BOSE-EINSTEIN CORRELATION FUNCTION 

The two-particle correlation function of two particles with 
four-momenta pi and p2 is given by the ratio of the two- 
particle number density, p2(/?i,p2), to the product of the two 
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single-particle number densities, pi{pi)pi{p2). Since we 
are here interested only in the correlation R2 due to Bose- 
Einstein interference, the product of single-particle densities 
is replaced by po{pi,p2), the two-particle density that would 
occur in the absence of Bose-Einstein correlations: 



R2{Pl,P2) 



PljPhPl) 

Poipl, Pi) 



(1) 



This P2 is corrected for detector acceptance and efficiency us- 
ing Monte Carlo events, to which a full detector simulation 
has been applied, on a bin-by-bin basis. An event mixing 
technique is used to construct po. This technique removes 
all correlations, not just Bose-Einstein. Hence, po is corrected 
for this using Monte Carlo lITlll . 

Since the mass of the two identical particles of the pair is 
fixed to the pion mass, the correlation function is defined in 
six-dimensional momentum space. Since Bose-Einstein cor- 
relations can be large only at small four-momentum differ- 
ence Q, they are often parametrized in this one-dimensional 
distance measure. There is no reason, however, to expect the 
hadron source to be spherically symmetric in jet fragmenta- 
tion. Recent investigations have, in fact, found an elongation 
of the source along the jet axis fH, [H, [13, H. While this 
effect is well established, the elongation is actually only about 
20%, which suggests that a parametrization in terms of the 
single variable Q, may be a good approximation. 

This is not the case in heavy-ion and hadron-hadron inter- 
actions, where BEC are found not to depend simply on Q, 
but on cornponents of the momentum difference separately 
JH [1^ I20I I21I I22I1 . However, in e+e^ annihilation at lower 
energy ll23ll it has been observed that Q is the appropriate 
variable. We checked this and confirm that this is indeed 
the case: We observe (Tl^ that R2 does not decrease when 
both — {pi — P2)^ and — {Ei —£'2)^ are large while 



q — qQ is small, but is maximal for Q — q 



0, 



independent of the individual values of q and ^o- The same 
is observed in a different decomposition: ~ Qt+Q^B^ 
where Qf = (pa — paf' is the component transverse to the 
thrust axis and Q^g = {p\\ ~ PxiY ~ (^1 ^£^2)^ combines the 
longitudinal momentum and energy differences. Again, R2 is 
maximal along the line Q = 0, as is shown in Fig.[T] This is 
observed both for two-jet and three-jet events. We conclude 
that a parametrization in terms of Q can be considered a good 
approximation for the purposes of this article. 
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FIG. 1 : R2 for two-jet events as function of the squares of the trans- 
verse momentum difference and the combination of longitudinal mo- 
mentum difference and energy difference. 



in. PARAMETRIZATIONS OF BEC 

With a few assumptions 10. S IIS], the two-particle corre- 
lation function, Eq. is related to the Fourier transformed 
source distribution: 



R2{pi,P2)=y[i+x\m)\^]{i+m 



(2) 



where f{x) is the (configuration space) density distribution of 
the source, and f{Q) is the Fourier transform (characteristic 
function) of f{x). The parameter y and the (1 +bQ) term 
have been introduced to parametrize possible long-range cor- 
relations not adequately accounted for in the reference sample, 
and the parameter X to account for several factors, such as the 
possible lack of complete incoherence of particle production 
and the presence of long-lived resonance decays if the particle 
emission consists of a small, resolvable core and a halo with 
experimentally unresolvable large length scales ll25il26ll . 
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FIG. 2: The Bose-Einstein correlation function R2 for two-jet events 
with the result of a fit of (a) the Gaussian and (b) the Edgeworth 
parametrizations, Eqs. ^ and respectively. The dashed line rep- 
resents the long-range part of the fit, i.e., "{{I +8Q). 



A. Gaussian distributed source 

The simplest assumption is that the source has a 
symmetric Gaussian distribution, in which case f{Q) = 



exp ijjQ - 



and 



R2{Q)=y[l+Xexp{-{RQf)]{l+5Q) 



(3) 



A fit of Eq. (|3]l to the data results in an unacceptably low 
confidence level. The fit is particularly bad at low Q values, 
as is shown in Fig.|2^ for two-jet events and in Fig.[3t for 
three-jet events, from which we conclude that the shape of the 
source deviates from a Gaussian. 

A model-independent way to study deviations from the 
Gaussian parametrization is to use jsl |23, IHl the Edgeworth 



expansion l'29'l about a Gaussian. Keeping only the first non- 
Gaussian term, we have 



R2{Q)=y{l+Xexp{-iRQf 



l + -H3{RQ)\)il+dQ) , 
(4) 

where K is the third-order cumulant moment and Hj,{RQ) = 
{y/lRQ)^ — 3V2RQ is the third-order Hermite polynomial. 
Note that the second-order cumulant coiTesponds to the ra- 
dius R. 

A fit of Eq. (|4|i to the two-jet data, shown in Fig.lJJ), is in- 
deed much better than the purely Gaussian fit. However, the 
confidence level is still marginal, and close inspection of the 
figure shows that the fit curve is systematically above the data 
in the region 0.6-1.2 GeV and that the data for Q > 1.5 GeV 
appear flatter than the curve, as is also the case for the purely 
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FIG. 3: The Bose-Einstein correlation function R2 for three-jet 
events with the result of a fit of (a) the Gaussian and (b) the Edge- 
worth parametrizations, Eqs. (O and {4}, respectively. The dashed 
line represents the long-range part of the fit, i.e., y{l + 62). 



Gaussian fit. Similar behavior is observed for three-jet events 
(Fig. |3j') and for all events. 



B. Levy distributed source 

The symmetric Levy stable distribution is characterized by 
three parameters: xo, R, and a. Its Fourier transform, f(Q), 
has the following form: 



f{Q) = exp iQ.xo 



(5) 



The index of stability, a, satisfies the inequality < a < 2. 
The case a = 2 corresponds to a Gaussian source distribution 
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FIG. 4: The Bose-Einstein correlation function R2 for two-jet events. 
The curve corresponds to the fit of the symmetric Levy parametriza- 
tion, Eq. l[6j. The dashed line represents the long-range part of the 
fit, i.e., y(1 + 8Q). The dot-dashed line represents a linear fit in the 
region Q > 1.5 GeV. 



with mean xq and standard deviation R. For more details, see, 
e.g., M- 

Then R2 has the following, relatively simple, form lf3lll : 



R2{Q)=y[l+lexp{-{RQr)]{l+dQ). 



(6) 



From the fit of Eq. (|6]l to the two-jet data, shown in Fig. |4] 
it is clear that the correlation function is far from Gaussian: 
a = L34±0.04. The confidence level, although improved 
compared to the fit of Eq. (O, is still unacceptably low, in 
fact worse than that for the Edgeworth parametrization. The 
same is true for three-jet events (Fig. |5]l and for all events. The 
values of a are 1 .39 ± 0.04 for three-jet and 1 .43 ± 0.03 for all 
events, respectively. 

Both the symmetric Levy parametrization and the Edge- 
worth parametrizations do a fair job of describing the re- 
gion Q < 0.6 GeV, but fail at higher Q. R2 in the region 
Q>1.5 GeV is nearly constant (« 1). However, in the region 
0.6-1.5 GeV R2 has a smaller value, dipping below unity |57], 
which is indicative of an anti-correlation. This is clearly seen 
in Figs. |4] and |5]by comparing the data in this region to an 
extrapolation of a linear fit, Eq. (|6]l with A- = 0, in the region 
Q> 1.5 GeV. The inability to describe this dip in R2 is the 
primary reason for the failure of both the Edgeworth and sym- 
metric Levy parametrizations. 



C. Time dependence of the source 

The parametrizations discussed so far, which have proved 
insufficient to describe the BEC, all assume a static source. 
The parameter R, representing the size of the source as seen 
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FIG. 5: The Bose-Einstein correlation function R2 for three-jet 
events. The curve corresponds to the fit of the symmetric Levy 
parametrization, Eq. The dashed line represents the long-range 
part of the fit, i.e., y(1 +8Q). The dot-dashed line represents a linear 
fit in the region g > 1.5 GeV. 



in the rest frame of the pion pair, is a constant. It has, how- 
ever, been observed that R depends on the transverse mass, 
y/E^ - pI, of the pions IH S. It has 
Il34i I35I1 that this dependence can be understood 



nit = 

been shown 

if the produced pions satisfy, approximately, the (generalized) 
Bjorken-Gottfried condition |l36l[37tl38ll3|,|40iySl], whereby 
the four-momentum of a produced particle and the space-time 
position at which it is produced are linearly related; 



(7) 



Such a correlation between space-time and momentum- 
energy is also a feature of the Lund string model as incor- 
porated in Jetset, which is very successful in describing de- 
tailed features of the hadronic final states of e+e^ annihila- 
tion. 

In the previous section we have seen that BEC depend, at 
least approximately, only on Q and not on its components sep- 
arately. This is a non-trivial result. For a hydrodynamical type 
of source, on the contrary, BEC decrease when any of the rela- 
tive momentum components is large llH Hit] . Further, we have 
seen that R2 in the region 0.6-1 .5 GeV dips below its values at 
higher Q. 

A model which predicts such Q-dependence while incor- 
porating the Bjorken-Gottfried condition is the so-called T- 
model, described below. 



1. The X-model 

A model of strongly correlated phase-space, known as the 
T-model ll42ll . explains the experimentally found invariant rel- 
ative momentum dependence of Bose-Einstein correlations in 



e+e^ reactions. This model also predicts a specific transverse 
mass dependence of R2, that we subject to an experimental 
test here. 

In this model, it is assumed that the average production 
point in the overall center-of-mass system, x = (t,rx,ry,r^), 
of particles with a given four-momentum k is given by 



In the case of two-jet events, 

d — t/oti , 



(8) 



(9) 



where nii is the transverse mass and z — yt — rz is the lon- 
gitudinal proper time [58]. For isotropically distributed par- 
ticle production, the transverse mass is replaced by the mass 
in Eq. while for the case of three-jet events the relation 
is more complicated. The second assumption is that the dis- 
tribution of x^'{k'') about its average, 5a(x^(^'') ~x^{k'')), is 
narrower than the proper-time distribution. Then the emission 
function of the T-model is 

S{x,k)= f dT:H(T:)5A(x-dk)pi(k) , (10) 
Jo 

where H{z) is the longitudinal proper-time distribution, the 
factor 8a {x — dk) describes the strength of the correlations be- 
tween coordinate space and momentum space variables and 
Pi{k) is the experimentally measurable single-particle spec- 
trum. 

The two-pion distribution, P2{ki,k2), is related to S{x,k), 
in the plane-wave approximation, by the Yano-Koonin for- 
mula lij: 



P2 



= j df'x\A'^X2S(x\,k\)S{x2,k2) 

•(l+cos([fci-fc2][xi-X2])) . (11) 



Approximating the function 5a by a Dirac delta function, the 
argument of the cosine becomes 



{kx-k2){xi -X2) - -Q.5{dx+d2)Q^ 



(12) 



Then the two-particle Bose-Einstein correlation function is 
approximated by 



^2(^1,^2) = 1 +XRei/^ ( ^ ) , 



(13) 



where //(co) = / dT//(T)exp(/a)T) is the Fourier transform of 
//(t). Thus an invariant relative momentum dependent BEC 
appears. Note that R2 depends not only on Q but also on the 
average transverse mass of the two pions, mt. 

Since there is no particle production before the onset of 
the collision, H{x) should be a one-sided distribution. We 
choose a one-sided Levy distribution, which has the charac- 
teristic function |31] (for a 7^ 1) 



//(co) = exp 



-^(At|co|)^ 



1 — / sign (co) tan 



(f)) 



-iCOTo 

(14) 
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where the parameter To is the proper time of the onset of parti- 
cle production and At is a measure of the width of the proper- 
time distribution. For the special case a = 1, see, e.g., IBOll . 
Using this characteristic function in Eq. (fTsl l yields 



1 



expl 



nil 



tan 



a7t\ 



2mt ^ 



2mt 



(1 + 86) 



(15) 



2. The x-modelfor average iiit 

Before proceeding to fits of Eq. ( fTSl l. we first consider a 
simplification of the equation obtained by assuming (a) that 
particle production starts immediately, i.e., Xo = 0, and (b) an 
average njt-dependence, which is implemented in an approx- 
imate way by defining an effective radius, R = \/Ai/ (2mt). 
This results in: 

^2(2) = y [1 + Xcos [(/Jag)'"] exp {-{RQf'')] (1 + 56) , 

(16) 



where is related to R by 



R 



2a 



/ajt\ 
:tan(-)7^ 



,2a 



(17) 



Fits of Eq. ( fTSI l are first performed with R^ as a free parameter. 
The fit results obtained, for two-jet, three-jet, and all events 
are listed in Table|T]and shown in Fig. |6]for two-jet events and 
in Fig. |7]for three-jet events. They have acceptable confidence 
levels, describing well the dip below unity in the 0.6-1.5 GeV 
region, as well as the low-Q peak. 

The fit parameters for the two-jet events satisfy Eq. ( fTTI l. 
However, those for three-jet and all events do not. We note 
that the values of the parameters a and R do not differ greatly 
between 2- and 3-jet samples, the most significant difference 
appearing to be nearly 3o for a. However, these parameters 
are rather highly correlated (in the fit for all events, the cor- 
relation coefficients are p{X,R) = 0.95, p(A-,a) = —0.67 and 
p{R,a) — —0.61, which makes the simple calculation of the 
statistical significance of differences in the parameters unreli- 
able. 

Fit results imposing Eq. (fTTl l are given in Table HIl For two- 
jet events, the values of the parameters are the same as in the 
fit with /?a free — only the uncertainties have changed. For 
three-jet and all events, the imposition of Eq. ( fTTT i results in 
values of a and R closer to those for two-jet events, but the 
confidence levels are very bad, a consequence of incompati- 
bility with Eq. (TU . an incompatibility that is not surprizing 
given that Eq. (|9]l is only valid for two-jet events. Therefore, 
we only consider two-jet events in the remaining sections of 
this article. 



3. The x-model with rrit dependence 

Fits of Eq. ( fTSI l to the two-jet data are performed in several 
nil intervals. The resulting fits are shown for several nii in- 
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FIG. 6: The Bose-Einstein correlation function R2 for two-jet events. 
The curve corresponds to the fit of the one-sided Levy parametriza- 
tion, Eq. ( 116b . The dashed line represents the long-range part of the 
fit, ;-.e., 7(1+8(2). 
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FIG. 7: The Bose-Einstein correlation function Ri for three-jet 
events. The curve corresponds to the fit of the one-sided Levy 
parametrization, Eq. ( 116b . The dashed line represents the long-range 
part of the fit, i.e., y(1 +§2). 



tervals in Fig. [8] and the values of the parameters obtained in 
the fits are listed in Fig.|9] The quality of the fits is seen to 
be statistically acceptable and the fitted values of the model 
parameters, a. To and At, are stable and within errors in- 
dependent of nil, confirming the expectation of the T-model. 
We conclude that the T-model with a one-sided Levy proper- 
time distribution describes the data with parameters To ~ fm, 
a ~ 0.38 ±0.05 and At « 3.5 ±0.6 fm. These values ai-e con- 
sistent with the fit of Eq. (fTSI l in the previous section, includ- 
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TABLE I: Results of fits of Eq. J16b for two-jet, three-jet, and all 
events. The uncertainties are only statistical. 



parameter 


2-jet 


3-jet 


all 


a 


0.42 ± 0.02 


0.35 ±0.01 


0.38 ± 0.01 


X 


0.67 ± 0.03 


0.84 ± 0.04 


0.73 ± 0.02 


R(fm) 


0.79 ± 0.04 


0.89 ± 0.03 


0.81 ±0.03 


(fm) 


0.59 ± 0.03 


0.88 ± 0.04 


0.81 ±0.02 


5 


0.003 ± 0.002 


-0.003 ± 0.002 0.003 ± 0.001 


Y 


0.979 ± 0.005 


1.001 ±0.005 


0.997 ± 0.003 


X^/DoF 


97/94 


102/94 


98/94 


confidence level 


40% 


27% 


37% 



TABLE 11: Results of fits of Eq. imposing Eq. il7b for two-jet, 
three-jet, and all events. The uncertainties are only statistical. 



parameter 


2-jet 


3-jet 


all 


a 


0.42 ± 0.01 


0.44 ± 0.01 


0.45 ± 0.01 


X 


0.67 ± 0.03 


0.77 ± 0.04 


0.69 ± 0.03 


R(fm) 


0.79 ± 0.03 


0.84 ± 0.04 


0.79 ± 0.03 


5 


0.003 ± 0.001 


0.010 ±0.001 


0.009 ± 0.001 


Y 


0.979 ± 0.005 


0.972 ±0.001 


0.973 ± 0.001 


X-^/DoF 


97/95 


174/95 


175/95 


confidence level 


42% 


10"'' 


10"^ 




0.5 1 
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FIG. 9: The fit parameters from fits of Eq. il5l to two-jet data for 
various intervals of mf 



IV. THE EMISSION FUNCTION OF TWO- JET EVENTS 



ing the value of R, which, combined with the average value of 
nil (0.563 GeV), corresponds to At = 3.5 fm. Just as in the fit 
of Eq. ( fTSI l. the parameters of the Levy distribution are highly 
correlated. Typical values of the correlation coefficients are 
P(A.,At) = 0.95, p(A.,a) = -0.67 and p(AT,a) = -0.9. 
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FIG. 8: The results of fits of Eq. l |15l ) to two-jet data for various 
intervals of nit. 



Within the framework of the x-model, we now reconstruct 
the space-time picture of the emitting process for two-jet 
events. The emission function in configuration space, S{x), 
is the proper time derivative of the integral over k of S{x,k), 
which in the T-model is given by Eq. ( fTOl i. Approximating 5^ 
by a Dirac delta function, we find 



S{x) 



A^n 
diXd?r 



mtr\ 

T / 



(18) 



To simplify the reconstruction of S{x) we assume that it can 
be factorized in the following way: 



S{r,z,t)=I{r)G{^)H{x) 



(19) 



where I{r) is the single-particle transverse distribution, G(r|) 
is the space-time rapidity distribution of particle production, 
and H{z) is the proper-time distribution. With the strongly 
correlated phase-space of the T-model, V[—y and r = pii/nii. 
Hence, 



G(ii) - Ny{y\) , 

3 

^ J Np^rnit/z) , 



Hr) . f^^^ 



(20) 
(21) 



where A', and A^^^ are the single-particle inclusive rapidity and 
pt distributions, respectively. The factorization of transverse 
and longitudinal distributions has been checked. The distri- 
bution of Pt is, to a good approximation, independent of the 
rapidity [11]. 

With these assumptions and using //(x) as obtained from 
the fit of Eq. (flST l (shown in Fig.fTOli together with the inclu- 
sive rapidity and pt distributions [11], the full emission func- 
tion is reconstructed. Its integral over the transverse distri- 
bution is plotted in Fig.[TT] It exhibits a "boomerang" shape 
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FIG. 1 1 : The temporal-longitudinal part of the source function nor- 
malized to the average number of pions per event. 



with a maximum at low t and z but with tails reaching out to 
very large values of t and z, a feature also observed in hadron- 
hadron |44] and heavy ion collisions ||45ll . 

The transverse part of the emission function is obtained by 
integrating over z and azimuthal angle. Figure [12] shows the 
transverse part of the emission function for various proper 
times. Particle production starts immediately, increases 
rapidly and decreases slowly. A ring-like structure, similar to 
the expanding, ring-like wave created by a pebble in a pond, 
is observed. These pictures together form a movie that ends in 
about 3.5 fm, making it the shortest movie ever made of a pro- 
cess in nature. An animated gif file covering the first 0.3 fm 
(10^24 gg^.) jg available Q. 



V. DISCUSSION 

BEC of all events as well as two- and three-jet events are 
observed to be well-described by a Levy parametrization in- 



FIG. 12: The transverse source function normalized to the average 
number of pions per event for various proper times. 



corporating strong correlations between configuration- and 
momentum-space. A Levy distribution arises naturally from 
a fractal, or from a random walk or anomalous diffusion [47], 
and the parton shower of the leading log approximation of 
QCD is a fractal H IM In this case, the Levy in- 

dex of stability is related to the strong coupling constant, as, 
byillll 



as = 



27t 



-a 



(22) 



Assuming (generalized) local parton hadron duality fsH 113, 
l55ll . one can expect that the distribution of hadrons retains the 
features of the gluon distribution. For the value of a found 
in fits of Eq. ( fT6] l we find as = 0.37 ± 0.04 for two-jet events. 
This is a reasonable value for a scale of 1-2 GeV, which is 
where the production of hadrons takes place. For comparison, 
from T decay, as(mx w 1.8 GeV) = 0.35 ±0.03 [56]. 

It is of particular interest to point out the mt dependence 
of the "width" of the source. In Eq. ( fTsT i the parameter as- 
sociated with the width is At. Note that it enters Eq. ( fTSl l as 
AxQ^/mt. In a Gaussian parametrization the radius R enters 
the parametrization as R^Qp'. Our observance that At is in- 
dependent of nix thus corresponds to R °^ 1 / \/Wi and can be 
interpreted as confirmation of the observance lf32ll33ll of such 
a dependence of the Gaussian radii in 2- and 3-dimensional 
analyses of Z decays. The lack of dependence of all the pa- 
rameters of Eq. ( fTsT i on the transverse mass is in accordance 
with the T-model. 

Using the BEC fit results and the T-model, the emission 
function of two-jet events is reconstructed. Particle produc- 
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tion begins immediately after collision, increases rapidly and pands outwards, which is similar to the picture in hadron- 
then decreases slowly, occuring predominantly close to the hadron interactions but unhke that of heavy ion collisions, 
light cone. In the transverse plane a ring-like structure ex- 
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